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1 Introduction 


The basis constructions considered in this paper involve elements from the 
theory of operator algebras. Since this may not be widely known to readers in 
approximation theory, we begin with a few remarks. 

The study of operator algebras breaks up in two parts: One the study of “the 
algebras themselves” as they emerge from the axioms, von Neumann algebras, 
and C'*-algebras. The other has a more applied slant: It involves “representa¬ 
tions” of the algebras. There is a close connection between the two parts of the 
theory: For example, representations of O^-algebras generate von Neumann al¬ 
gebras. It was realized in the last ten years (see e.g., [BJP96], [DuJoOS]) that 
certain families of representations are useful in basis constructions in harmonic 
analysis, in approximations, in signal/image analysis, and more generally in 
computational mathematics. The bases in question may typically be built up 
from representations of an especially important family of simple C^-algebras, 
known as the Cuntz algebras [Cu77]. These Cuntz algebras (see Lemma 2.2 
below) are denoted 02 , 03 ,including Ooo- 

The connection to Cuntz algebras Ojv is further relevant to the kind of dynami¬ 
cal systems built on iterated branching-laws, with the case of Ojv representing 
A^-fold branching. The reason for this is that if N is fixed, Ojv includes in 
its dehnition an iterated branching, taking the form of subdivision, but now 
within the context of Hilbert space; so we generate subdivisions into orthogo¬ 
nal families of subspaces of the initial Hilbert space. In this paper, we follow 
up on a certain probabilistic aspect of this construction. 

The OjY point of view is especially well suited to basis constructions in such 
contexts as wavelets and fractals since they naturally involve the same kind of 
sub-division. Our starting point is an initial Hilbert space H, where H may be 
or Ti may be for some fractal measure p; see Section 2 below. 

The more successful bases in Hilbert space are the orthonormal bases ONBs, 
and we shall consider a certain computational algorithm for generating them. 

A further reason the subdivision schemes in Hilbert space are useful is that 
the more familiar Fourier wave functions are periodic, and so not localized. 
Moreover these existing Fourier tools are typically not friendly to algorithmic 
computations. By a local (See Dehnition 1.1) construction we mean an algo¬ 
rithm which begins with a hnite family of functions (often one or two) having 
a hxed compact (i.e., local) support, and a procedure allowing assigned scaling 
and translation operations. As is popular in the sub-band approach to wavelets 
and wavelet packets in ID, the scaling is typically in powers of a hxed base, 
i.e., it could consist of powers where N is hxed {N > 1) and where j varies 
over Z, stretching and squeezing the support. 
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The main result in this paper concerns bases in the Hilbert spaces p) de- 

hned from measures p arising as equilibrium measures (also called self-similar 
measures) for iterated function systems (IFS). However, our results are spelled 
out in more detail for T^((0,1); Lebesgue), where the classical Walsh system is 
a special case. Our construction uses ideas from dynamical systems and oper¬ 
ator algebras (specihcally representations of the Cuntz algebras). The Cuntz 
algebra [Cu77] is used in the construction of bases for geometric structures 
with self-similarity such as iterated function systems (IFS). Recall that for 
a hxed hnite N, the Cuntz algebra Oat is generated abstractly by N isome¬ 
tries. In representations of Otv on a concrete Hilbert space Ti, the resulting 
isometries Si, say, have orthogonal ranges which form a partition of unity in 
the particular Hilbert space which carries the representation, in the sense 
that the identity operator is written as a sum of the N projections SiS* 
onto the respective ranges SiH. Since the subdivision process can be iterated, 
this idea has already proved useful in understanding orthogonal families in 
Hilbert spaces built on IFSs; see, e.g., [DuJo06b,JoPe96,JoPe98]. Our analy¬ 
sis here uses such particular On representations in combination with certain 
graph-theoretic considerations. In addition to the reference [Cu77], the Cuntz 
algebras On and their representations are reviewed in [Jor06, Sections 7.6 and 
7.7]. The book [Jor06] also includes additional motivation, and [Jor06, Chap¬ 
ter 4, p. 69, and Section 9.4] cover details for an IFS-family of Cantor systems 
and their self-similar measures. 

Applications to Brownian motions are intended but postponed to a later paper. 
The connection between Brownian motion, Cuntz algebras, and IFSs is treated 
in the literature, for example in [Jor06, pp. 56-57, 151, and 203]. 

There has been a recent increased interest in basis constructions outside the 
traditional setting of harmonic analysis. The setting which so far has proved 
more amenable to an explicit analysis with basis functions involves a mix of 
analysis and dynamics, and it typically goes beyond the standard and more fa¬ 
miliar setting of orthonormal bases consisting of Fourier frequencies. The con¬ 
text of frames in Hilbert space (see, e.g., [ALTW04], [BoPa05], and [BPS03]) 
is a case in point. 

As is well known, the classical setting of Fourier analysis presupposes a choice 
of Fourier frequencies, or Fourier trigonometric basis functions. However, this 
unduly limits our choices, and the applications: As is well known, Fourier’s 
basis functions are less localized, and the computational formulas typically 
are not recursive. There are now alternative dynamical approaches which are 
recursive, and at the same time are amenable to harmonic basis constructions; 
see, e.g., [Dut05], [JoPe98], [DuJo05], and [DuJo06a]. Moreover, these recursive 
models arise in applications exhibiting a suitable scale-similarity. Their con¬ 
sideration combines classical ideas from inhnite convolution with ideas from 
dynamics of a more recent vintage. 
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In this paper, we introduce a recursive scheme for a basis construction in 
the Hilbert space -h^(0,1) which is analogous to that of Haar [HaalO] and 
Walsh [Wal23,Chr55]. While computationally efficient, these more traditional 
approaches limit the choices of functions too much, often to step functions, or 
at any rate to functions that have a limited number of derivatives. 

We begin here with a certain dual system of axioms for reflection symmetries 
for functions on the unit interval X = [0,1]. We then show that up to this 
reflection symmetry, we get recursive algorithms which in turn produce inhnite 
families of orthonormal bases localized in X^(0,1). Our scheme is as versatile 
as the more traditional spline constructions; see, for example, [Wic94]. But 
while singly generated spline systems do not produce orthonormal bases, each 
of our present algorithms does. Moreover our scheme is adapted to the hxed 
unit interval X = [0,1] while the more traditional wavelet-based wavepackets 
are designed for the construction of orthonormal bases in X^(M); see [CoWi92]. 
And if the starting functions are of compact support, the size of the support 
reaches outside the unit interval [0,1]. 

The more traditional approaches to basis algorithms further have limited the 
libraries of functions to be used at the initial step of the recursion. We get 
around this here by identifying a set of symmetry conditions that may be 
imposed on two functions -ipo and in the Hilbert space X^(0,1); see Fig. 
1 for an illustration in the simplest case. Our algorithm is then based on a 
certain matrix scaling and subdivision applied to these two functions. Hence 
our starting point is different from the more traditional one which begins 
with a scaling identity, masking coefficients, and a so-called father function (fo 
which solves the corresponding scaling identity; see [Dau92]. 

Our justihcation for the term “wavelet” in connection with the present basis 
is threefold: 

(a) Our functions are localized in a sense which will be made clear. 

(b) Our construction is recursive. 

(c) Our algorithm for constructing orthonormal bases starts with a prescribed 
and carefully selected hnite system S of functions. Two operations are applied 
recursively to S, scaling and reflection. But note that as the algorithm runs, 
the reflections are scaled as well. 

Definitions 1.1. A family of functions (typically an ortho normal basis) {ipn '■ 
n > 1} in X^(0,1) is said to be local if for all e > 0 there exists an N{e) > 1 so 
that the closed linear space spanned by {tjjn : I < n < iV(e)} admits a hnite 
subfamily of orthogonal functions {(pn : I < n < M(e)} total for the subspace 
spanned by {-ipn : I < n < iV(e)} with the restriction that the Lebesgue 
measure of the support of each ipn is less then e. 
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Fig. 1. The first thirty-two functions in the sequence (pn 

To understand the two discrete operations which underlie our construction, 
it is helpful to review a fundamental feature of wavelets in the non-standard 
setting of iterated function systems (IFS); see also [JoPe94,JoPe96] and Sec¬ 
tion 4 below for additional details. In the simplest of settings, i.e., that of the 
unit interval, the two discrete operations going into our algorithms are that of 
iterated scaling by 2, i.e., the system x 2'^x, modulo 1, for m = 0,1, 2,..., 
and mid-point reflections. Both operations have analogues for more general 
IFSs, and we discuss these generalizations below. 
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Fig. 1 (cont’d). 


Our use of the terms “reflection” and “reflection symmetry” is related to, but 
different from, the one studied in, for example, [Jo0198]. The main difference 
is that in our present context, the Hilbert-space inner product is preserved 
after the reflection, while in [Jo0198] it is changed, i.e., it is subjected to a 
certain renormalization. 
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2 Iterated function systems 


To understand our results it is useful to consider a slightly more general setup: 
Let d be a fixed dimension, and consider a given finite set S of affine and 
contractive mappings r* : There are several interesting limits in 

the literature, arising from iteration of such a system S = (tj). The accepted 
terminology is “affine Iterated Function System (IFS)”. See [HutSl], [Jor06], 
and formulas (2.1)-(2.2) below. 

Following Cantor’s middle-third construction, note that one limiting object 
derived from S results from recursive iterations of the individual maps in S'; 
it is an attractor which takes the form of a compact subset X{= X{S)) of 
The set X{S) often has fractal like properties: And in one dimension {d = 1) 
the deleted-middle-third Cantor set is an example, but the unit interval is 
one too. The other iteration limits in this context take place in the family of 
probability measures on M'’*. However when S is given, the limit measure in 
question (called equilibrium measure) depends further on a chosen assignment 
of probability weights: It turns out that, given S, and given a fixed assignment 
of probabilities (pj) to the TjS, there is a unique equilibrium measure p = fis,p- 
If Pi > 0 for all i, then the support of the measure fis,p is X{S). Even in ID 
when X [S) may be the unit interval, there is a variety of measures arising from 
the second limit construction other than the restricted Lebesgue measure. 

In this paper, we are interested in localized orthonormal bases (ONBs) in 
L‘^{^s,p)- If turns out that our main issues may be best presented in the case 
of d = 1, and in the special case where the weights are uniform, i.e., p* = 1/N 
where N is the number of maps Ti from the initial system S. In fact, as noted 
in [DuJoObb], the Hilbert space L‘^{fis,p) does not have ONBs of complex 
exponential Fourier bases in the case of non-uniform weights. 

While the gist of our paper is for the unit interval, we wish to add that our 
construction in fact works in a more general context, that of iterated function 
systems (IFS) [Hut81,JoPe94]. However, it is easier to get an overview of the 
totality of admissible bases in the special case of the unit interval X = [0,1]. 

We consider decomposition theory for the Hilbert space of square-integrable 
functions on the unit interval, both with respect to Lebesgue measure, and 
also with respect to a wider class of self-similar measures p [JoPe96]. That 
is, we consider recursive and orthogonal decompositions for the Hilbert space 
X^(p) where p is some self-similar measure on [0,1]. 

We now turn to the technical details needed in our discussion of IFSs in 
Sections 3 and 4 below. For the axiomatics of IFSs, see, e.g., [HutSl] and 
[JorOSj. A finite iterated function system is determined by a finite system of 
contractive endomorphisms r = (r^) in a compact metric space X. We shall 
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consider here a finite system consisting of N endomorphisms of a compact 
space X. As we will see, one of the endomorphisms will be singled out, and 
it will be convenient to index it by zero, i.e., the first map is Tq. When such 
a system r is given, it is then known that there is a unique Borel probability 
measure n on X such that 


/i 



o Ti 


-1 


( 2 . 1 ) 


or equivalently, 

fix) dfiix) = J /(A(a^)) dnix) for all / e C(X). 

The measure n is called the Hutchinson measure, and it is also called the 
balanced invariant measure of r. In the case when X is [0,1], the unit interval, 
then the two maps tq and ri may be taken to be a; —a;/2 and a; —> (x + l)/2, 
respectively, and fi is then the standard normalized Lebesgue measure on [0,1]. 
For general IFSs (X, r), the starting point for our ONB construction in the 
Hilbert space L‘^{X,fi) is then a specified finite set of functions [ipi) which 
satisfy a certain reflection symmetry which we proceed to describe in detail in 
Sections 3 and 4 below. 


Let {X,d) be a compact metric space. Let a: X — X be an endomorphism 
such that the number of elements in a~^{x} is equal to N for all x E X, where 
X, 2 < X < cx), is fixed. Iteration of branches 

W^({a;}) = {y eX \ aiy) = x} 

then gives rise to a combinatorial tree. If 

^ = (i^i, 1^2, • • •) G H = {0,1,..., X — 1}^, 


an associated path may be thought of as an infinite extension of finite walks 

■ ■ ■ X2X12) 

with starting point x, where (rj), i = 0,1,..., X — 1, is a system of Borel 
measurable inverses of a, i.e.. 


a o n = idx, 0 < i < X — 1. 

We assume that each (rj) is contractive, i.e., there exists a constant 0 < c < 1 
such that 


diTiX,Tiy) < cdix,y), 


x,y eX, i e {0,1,... ,X - 1}. 



This contractive condition ensures that there exists a Borel probability mea¬ 
sure p on X such that 

^ N-l 

= X E 

i=0 

Before we state our hrst simple lemma we also set 

ai = a on rj(X). 

For each j G {0,1,..., iV — 1} we dehne a linear operator (Sj) by 

( 2 . 2 ) 

k=0 

for aA\ f en = 

Lemma 2.1. If the operators Sj in L^{X,fj,) are given as in (2.2), then the 
formula for the corresponding adjoint operators S* is as follows: 

1 -k 

•s;/ = K7 E ° n (2.3) 

k=0 


PROOF. We leave the details to the reader. They are based on the Tj- 
equipartition property for the measure p. See (2.1), and also [JoPe94]. □ 


In the next lemma we show that the operators Sj from (2.2) dehne a represen¬ 
tation of the Cuntz C*-algebra. As noted in [JoPe94,JoPe96], this C*-algebra 
has found a variety of uses in iterated function systems (IFS), and in approx¬ 
imation theory. Its use in analysis and in physics was initiated in the paper 
[Cu77]. For each IFS with N endomorphisms, there is an associated represen¬ 
tation of the Cuntz algebra with N generators Sj. These generators are 
also called the fundamental isometries. (When we say “isometry”, we are here 
referring to the Hilbert space L‘^{X,fi).) From [Cu77] we further know that 
specifying a representation of On is equivalent to specifying a system of N 
fundamental isometries. 

Lemma 2.2. The system of operators {Sj) from (2.2) defines a representation 
of the Cuntz algebra On- We write {Sj) G Rep((P(i, Ti), he., 

N-l 

S*Sk = E SkS: = (2.4) 

fc =0 


PROOF. Again, we leave the details to the reader. The argument uses the 
previous lemma combined with the axioms for the system (tj) and the asso¬ 
ciated measure p, outlined above. 
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Returning to formula (2.2) and setting 


(x) = Y. (a:), 

k=0 

we see that (2.2) may be rewritten in the form 

{Sjf) {x) = rrij {x) f (a (x)) for / G (X, /i), x E X. (2.5) 

Hence using an idea from [BJMP05], we then note that the Cuntz relations 
(2.4) for the operators {Sj) in (2.5) are equivalent to the assertion that the 
N X N matrix function U given by 

\ N-l 

ruj (Tfc (x)) 

/ j,k=0 

takes values in the group Utv (C) of all iV x iV unitary matrices. 

We now check that UU* = Iw- Substituting mj{Tk{x)) = into the 

formula for the matrix product, we get 

_ 1 ^-1 

E u,,u {x) Ui,k {x) = ^Y = hi ■ □ 

fc =0 fc =0 

Remark 2.3. There is a variety of representations of the Cuntz algebra Om-, 
other than those that come naturally from IFSs (as in (2.2) and the lemma). In 
fact by a theorem of Glimm [Gh60], the set of equivalence classes of irreducible 
representations oi On cannot be parametrized by a Borel cross section. 

Our focus in this paper is the search for orthonormal bases (ONBs) in the 
Hilbert space T^(/i) associated naturally with a hxed contractive IFS. Once an 
ONB is chosen, it may be used in the analogue-to-digital (A-to-D) conversion 
of signals. Conversely, for a given signal (read, / G L^(/i)) there is a variety 
of ONBs, and a choice must be made; for example, we may wish to minimize 
the entropy (information loss), 

[f) = f)\^\n\{'tl)^,\ f)h 

k 

where {'ipk) is the chosen ONB. (Here, we shall normalize /, i.e., assume 

ll/llL2(/i) = !•) 

We make use of representations in two ways: (1) as a hrst step toward the 
determination of an ONB; and (2) for entropy computations, even when an 
ONB is not known. 


c/(U = w>w)k:.k=;^( 
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The details on (1) will be presented in the next section; and we now briefly 
discuss (2). 


As illustrated in the lemma below, our use of the fundamental isometries {Si) 
in a particular representation of in some Hilbert space Ti. immediately 
yields scales of mutually orthogonal subspaces of H, or equivalently, mutually 
orthogonal projections which constitute partition of the identity operator II 
in Ti. Example: the projections P* := SiS* are orthogonal and satisfy 

As stressed in for example [KriOS] and [JKK05], the use of subspaces as op¬ 
posed to vectors is significant for the design of quantum algorithms such as 
quantum error-correction codes. The reason is that direct operations on in¬ 
dividual vectors in the Hilbert space H typically would destroy the quantum 
states on which the algorithm is operating. The quantum states (or qubits) 
might represent polarized photons. (Recall that in the conventions of quan¬ 
tum theory, unit-vectors in the underlying Hilbert space represent quantum 
states!) Our use of representations of lets the algorithms act on subspaces 
in Ti, and as a result leave the individual quantum states intact. The intrinsic 
orthogonality of the subspaces is what yields quantum channels; see [JKK05]. 


Lemma 2.4. Let {Si) G Rep((9 at, Ti) be a representation of On acting on a 
Hilbert space Ti. Let M.k denote the set of all multi-indices J = (ji j 2 • • • jk), 
where each ji is m {0,1,..., iV — 1}. Then for each k, 

Sj := Sj,---Sj^, JeMk, 

is a representation of On^; in particular, Pj := SjS} yields a commuting 
family of mutually orthogonal projections in TL, i.e., 

PjPj' = dgj'Pj, J G Mk, 

and 

E Pj = ^n- 

J&Mk 

For f in TL, ||/|| = 1, the entropy number 

JeMk 

satisfies 

^k+l if) = ^1 (/) + E ) 

where Pi = SiS* and 

N-l 

O (/) :=-EllE/ll'lnil^;/lE 

i=0 
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PROOF. The steps in the proof of the lemma are straightforward, and we 
leave them to the reader. □ 


Example 2.5. Let X = [0,1] and let a: X ^ X he the endomorphism 
a{x) = 2x mod 1. Here To{x) = | and Ti{x) = are two maps for which 
a o Ti{x) = X, i = 0,1. Let ao{x) = 2a;X[o,i] and ai{x) = {2x — l)x[i 4 ]; then 


Sof = f oao + f oai, 
Sif = f oao - f oai, 


Sy = °rQ + f oTi), 
SU = ^(/°T-o-/ori). 


We choose unit vectors if, tjj so that Slf = 0 and S'q'^ = 0. By Cuntz’s relation 
(2.4) we get S^S^f = f and SiSlip = Thus we get {f \ ip) = h. hiy our 
construction. 

Slip = 0 


if 


i.e.. 


1pOTo = -IpOTi, 
'X + 1 


V- (I) = -i> p 


or equivalently. 


iP (x) = -Ip =-ip[x + 


for all 0 < a; < -. 
- - 2 


Lemma 2.6. Let f be an element o/L^(0,1). Then the following statements 
are equivalent: 

(a) Slf = 0; 

(b) f{x) = -/ Q + for allt)<x< 

(c) f{x) = Y.(^n cos(27r(2n + l)a;) + bn sin(27r(2n + l)a;). 

n>0 

Moreover SiSl is the projection onto the closed subspace 

|/gL 2(0, 1) /(x) =-/Q+ a;^ , 0 < a; < ^ I . 


Similarly the following statements are equivalent: 


(a) Slf = 0; 

(b) f{x) = + 


for alio < X < 


1 

2 ’ 


(c) fix) = E cos( 27 r( 2 n)a;) + 6 nSin( 27 r( 2 n)a;). 


n>0 
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Moreover SqSq is the projection onto the closed subspace 

I / G L2(0, 1) f{x) = /Q + a;^, 0<a;<^|. 


PROOF. Equivalence of statements (a) and (b) is obvious. That (c) implies 
(b) is routine as 

cos(a; + (2m + l)7r) = — cos(a;) 

and 

sin(a; + (2m + l)7r) = — sin(a;) 

for any m > 0. We will prove now that (a) implies (c). To that end, for any 
n > 0 we set 

Cnif) = / cos{2TTnx) f (x) dx 
Jo 

and 

Snif) = / sm{2TTnx)f (x) dx. 

Jo 

A simple computation shows that c„(S'q/) = C 2 n(/) and SniS^f) = S 2 n{f) for 
all n > 0. Thus C 2 n(/) = S 2 n{f) = 0 if 5 * 0 / = 0. The last statement follows 
since by the Cuntz relations, S^f = 0 if and only if SiS^f = f. 

For equivalence of statements in the second set, we note that by the Cuntz 
relations, { f \ g) = 0 whenever S'q/ = 0 and S^g = 0. Thus, the equivalence of 
the second set of statements follows from that of the hrst set of statements. □ 

Notation 2.7. We shall need the subspaces 

}Co:= {'ipl = 0} 

and 

{Sr^p I S(/) < 1}. 

The symbol S(/) in the second subspace { | S(/) < 1}, S(/) = is 

dehned as all multi-indices I such that S(/) < 1. It is motivated as follows: 
In building ONBs, the aim is to start with a conveniently chosen function -ip, 
and then to construct the rest from recursively applying monomials in the 
generators Si (chosen from a particular On representation). This notation 
allows us to keep track of the combined system of relations in step-size of 
length one. 

For any ip with Sp^tp = 0, we have 


%p{x)'ip{p2x — 1) dx = 


I (^2 W 

1 

- 'ip{x)'ip{2x) dx. 
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Thus we have 


{ip \ Slip) = / ip{x)ip{2x) dx — ip{x)ip{2x — 1) dx 


= 2 


ip{x)ip{2x) dx. 


More generally for any two elements ip^^\ ip^‘^'^ of the subspace 


we have 


}Co:= {ip \ SqiP = 0}, 



2 ip‘^^\x)ip^‘^\2x) dx. 

Jo 


Any vector ip E H such that S^ip = 0 and {ip \ SiSp^ip) = 0 for all m > 0 
is called a generating vector for the closed linear span of the vectors { Sjip \ 
S(/) < 1} (see Notation 2.7 above and Lemma 2.10 below for the notation 
used here). A subspace /C of /Co is called a basis space for generating vectors if 
it is a maximal family of vectors that satishes the following mutnal relation; 

{ip \ SiSlP^iP') = 0, m>0, 

for all Ip, Ip' G JC. Existence of a snch a maximal family of vectors follows by 
Zorn’s lemma. It is simple to note that /C is a subspace of /Cq. 

Definition 2.8. A maximal subspace /C as above will be called a basis space. 

We now turn to a concrete representation of the subspaces /Co and /C in L^(0,1) 
which were ontlined above. 


Our identihcation of subspaces /Co and associated orthonormal bases is by a 
certain algorithmic procedure. Below we illustrate our recursive construction 
in one particular example: In our construction we begin with the representa¬ 
tion of O 2 from Lemma 2.2, and we give a natural and orthogonal subspace 
decomposition of the Hilbert space L^(0,1). Using this, we then show how 
an associated recursive basis may be realized. In our example, we start with 
a family of sine fnnctions, normalized to have period one. We then aim for 
an orthonormal basis (ONB) when the Hilbert space L^(0,1) is dehned from 
the restriction of Lebesgue measure to the unit interval X = [0,1]. (Other 
self-similar measures will be considered later!) Our example will further serve 
to illustrate the reflection operations which we will encounter later in a more 
general context of self-similar systems. 

As outlined before, the idea is to start our recursion from two prescribed 
fnnctions ipo and ipi. Here we take ipo to be the constant function “one” on 
X = [0,1]; and we choose ipi{x) := s{x) := sin(27ra;). The recursion will be 
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as outlined above: The idea is to recursively determine the pair of functions 
t/’ 2 n+i) from Sampled subdivisions of 'ipn for each n > 1. (Note that our 
recursion does not begin with n = 0.) Notation: Set Sn{x) := sin(27rna;), 
n = 1,2,.... 

Lemma 2.9. For every odd integer n, the function Sn is in JCq. Moreover, if 
n is even, then S^Sn is non-zero in L^(0,1). 


PROOF. We begin by setting s{x) = sin(27ra;), and for any integer n > 1 we 
also set Sn{x) = s{nx) = sin(27rna;), a; G [0,1]. As 


'1 + a;\ fl-\r X 

Si I —:— I = s 


^ = sin(7r(l + x)) 


X 

= -sin^TTXj = -Si ( - ) , 


we have si G /Co- For any odd integer, i.e., n = 2m + 1, we check that 

—j = sin(7r(2m + 1)(1 + x)) = sin(7r(2m + l)a; + (2m + l)7r) 
= — sin(7r(2m + l)a;) = —s„ ^ . □ 


We shall need the following additional facts about the functions For any 
two integers m, n > 0, we have 

1 

( Sm I 5*15^ ) 2 / Sm{x')Sn{2x') dx 

Jo 
1 

= / [cos(27r(m — 2n)x) — cos(27r(m + 2n)x)] dx 
Jo 


sin(27r(m — 2n)x) sin(27r(m + 2n)x) 


27r(m — 2n) 


27r(m + 2n) 


= 0 


for m 7 ^ 2n. For m = 2n, we also check that the integral is 0. Since = 
S 2 ^n, we also get {sm I SiS^Sn) = 0. Hence, {<S 2 n+i(a^) | n > 0 } is a family 
of orthonormal vectors in a basis space. One natural question that we face 
now: Is it a maximal family, i.e., is it a basis space? We answer this in the 
affirmative in the remaining part of this section. 

Lemma 2.10. Let ip be a unit vector such that 

SqiP = 0 and {ip \ S'i(S'o)™'t/’) = 0 for all m > 0, 

but 

%p I ^i(^o)”''^i(5'o)’”^) 7 ^ 0 for allm,m' ^ 0. 
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Then it follows that the family of vectors 




where 


and 


k 


I = {iii2 ... 4), ik e { 0 , 1 }, 


is a maximal family of orthonormal vectors. The closed subspace TLipp) spanned 
by the vectors {Sif: \ |J| < cx)} is invariant under both of the operators Sq 
and Sq. 


PROOF. Since Sqi/j = 0, it is simple to verify by the Cuntz relations that 
(S'o )™'4 is orthogonal to (S'o )”4 for m 7^ n > 0 , where by convention (S'o)° = 1. 
As S'J'S'o = 0 and = 0, we also check that Si{So)^ip is orthogonal 

to (S'o)”'S'i(S'o)”4 and (S'o)”V for all m' 7 ^ n' and n,n' > 0 . Thus we are left 
to check for m' = n'. In such case orthogonality follows by our hypothesis that 
(4 I 5'i(S'o)™'4) = 0 for all m > 0 . It is clear that the vector space generated 
by these vectors is both Sq- and S'g-invariant. The maximal property is also 
evident. See Section 3 below for details. □ 


Terminology: M. will denote the set of all hnite multi-indices. 

Proposition 2.11. Let {Si) be the irreducible representation of O 2 as in Ex¬ 
ample 2.5, and let 1C be a basis space. For a vector ip use the notation ^{{ip) 
as in Lemma 2.10 for the closed subspace spanned by the vectors { Sjip \ |/| < 
00 }. Then the following hold: 

(a) For each unit vector ip & 1C, the vectors in the family { Sjip \ S(/) < 1} 
are orthonormal, i.e., norm one, and mutually orthogonal. 

(b) For any two orthogonal unit vectors ip^S^^C^) ^ TL{ip^^'^) is orthogonal 

to TL{ip^‘^^). 


PROOF. Proof is routine as in Lemma 2.10. □ 

Lemma 2.12. {Haar-Walsh, Fig. 1 above, [HaalO,Wal23]) Let Ti = L^(0,1) 
and let ipo = X[o,i]- Then the recursive system 

<~P2n{x) = iPn{‘2x) -f iPni^X - 1), 

T 2 n+l{x) = (Pn{‘2x) - (Pn{2x - 1 ) 

for n > 0 defines an orthonormal basis for TL. 
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PROOF. We set ^o, S'! as in Example 2.5. We consider the Hardy space 
given by 


/(^) = i: 

n>0 

for {Cn) e P and ||/||^ = 


c tP' 

z eT^ = [z eC\ 

kl = 1}, 

X]n>0 

c„p. We also set en{z) 

= z^ and 


Sof{z) = f{z^), 

Sif{z) = zf{z^) 

for all / G 7-f+ and z eT^. We have S^^po = po and ^qCo = cq. Now, dehning 
W: 7-f+ —> Td by Wen = v^n, we verify that 

WSi = SiW for i e {0,1}. 

If n = ji + ^22 + ■ ■ ■ + Jfc2^ where jr E {0,1} is the dyadic representation of 
an integer n > 0, it follows that Wcn = W■ ■ ■ Sj^^eo = Sj^Sj^ ■ ■ ■ Sj^po = 
Pn, and the result follows. In the last step we are using the fact (details in 
[Jor05]) that the two O 2 representations are irreducible; so the intertwining 
operator W is a constant times a unitary. □ 

Definitions 2.13. (a) For p E'H \ {0}, set 

K (pp) := min{ K E M. \ s.t. ( Sj^ip \ ) = 0 V J* ^ ii', Ji 7 ^ J 2 , 

and 3 J ^ K s.t. ( Sjip \ SkP ) 7 ^ 0 }. 

(b) Further, define 

7i (pp) := closed span { SppSjpj \ m E Nq, J ^ K (pp) } . (2.6) 


(c) A family of vectors (ppn)i \\Pn\\ = 1, is said to be maximal and orthogonal if 
(i) the corresponding subspaces Ti (ppn) are mutually orthogonal, (ii) Ti ippn) -L 
p, and (hi) they are not part of a bigger such family. 

Theorem 2.14. Let 1C be a basis space, and let {pn \ n > 1} be an orthonor¬ 
mal basis for JC. Further, let Tin be a maximal family of subspaces associated 
with {pn) in the sense of Definition 2.13(c). Then 

n = Cp^nn. 

n>l 


Before turning to the proof, we will need some preliminaries which are included 
in the remarks below. The proof will then be resumed at the end of Section 3 
below. 
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The reasoning in the proof of the theorem has the following two parts in rough 
outline: Firstly, the argument for orthogonality of the spaces and the vectors 
which go into our basis construction is largely combinatorial, and it is sketched 
in Remarks 2.15 below. 

Secondly, we must prove that the functions which are produced by the al¬ 
gorithm form a total family in T^(0,1), i.e., that these vectors span a dense 
subspace in T^(0,1). Recall our algorithm starts with the two functions 'ipQ = 
ip = the constant function “one”, and = si(a;) = sin(27ra;). Using the 

sine functions from Lemma 2.9, and our representation of O 2 from Lemmas 
2.2 and 2.10, we then organize a system of orthogonal functions in each of 
the closed subspaces from the conclusion of the theorem. Our assertion 
is that this orthonormal family of vectors is total. Our argument going into 
the proof of this is structured as follows (details in Section 3): Suppose some 
/ in L^(0,1) is in the orthogonal complement of the family. We then show 
that both vectors / and S^f must have Fourier expansions consisting only 
of cosine functions. Because of the reflection built into the operator S^, we 
conclude that this is only possible if / is zero. 

Remarks 2.15. Our algorithm may be applied both to existing wavelets, and 
to new ones as well. Consider for example the Haar-Walsh sequence 

{iPn\n>0}CL^{0,l) 

dehned as in Lemma 2.12. For any multi-indices J = {jij 2 ■■■ jn) ^ M., 
ji G {0,1}, consider the orthogonal family Sjip and set 

Mev = {J e M \ S( J) = 0, 2,4, 6 ,... } 

and 

Ml = { Jg Af I S(J) < 1}. (2.7) 

Here a number N is fixed and our multi-indices are built from the alphabet 
{0,1,..., iV — 1}. Further we dehne S(J) to be ji- 

For any J E A4ev, set 

Tii^Sjip) = closed span{ Sxj'ip \ K G Mi }. 

We claim that 

®Cipo = L2(0, 1) (2.8) 

J 

and the vectors appearing in (2.8) are orthonormal. Note that the first few 
terms in the system of closed subspaces in (2.8) are 

'H{ipi) © 'H{ipj) © 7i{ipii) © 'H{ipiz) © • • • . 
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Specifically, 


if 7 = Sli), 
fii = sISq^p, 
fi^ = SiSqSiiP, 


We now turn to our basis constructions. 

In brief outline (Proof of Theorem 2.14 Part 1): The starting point is a hxed 
IPS with N endomorphisms, and an associated representation {Si) of (Pat in 
L^{f). Although it would seem natural to begin with cyclic subspaces in T^(/i), 
care must be exercised in order to guarantee orthogonality. Our starting point 
will be two orthogonal (and carefully selected) vectors f and ip satisfying 
S^f = f, and SqiP = 0. (Note that orthogonality of f and ip is implied by 
these relations.) Then using the representation, new vectors f are constructed 
to make up part of an ONB. The new vectors result from applying operator 
monomials in the generators Si to ip. So operator monomials are applied to the 
single vector ip, and the monomials are selected with view to orthogonality. 
The various new orthogonal vectors f are assigned subscripts according to 
the particular operator monomial which is applied to ip. The process is then 
repeated inductively on additional vectors ip as needed for creating a full ONB. 


3 Irreducible representations of On 


Let W G N, W > 2, and let Ti be a separable Hilbert space. Let {Si\0<i<N} 
be a system of isometries in Ti which dehne an irreducible representation of 
the Cuntz algebra On- We shall further assume that there is some f such 
that II99II = 1 and SQf = f. (Note that then Aq 99 = 99 as well.) 

Set M. = M-n (= the set of all hnite multi-indices J = (jo ji • • • ip) where 
ji G {0,1,..., W — 1}). Since 

(jo ji • • • ip) I—^ jo + iiN + ■ ■ ■ + ipN^ 

dehnes a bijection oiM .n onto No = {0, 1 ,... }, j^Tv acquires an order induced 
from the natural order on Nq. 

Theorem 3.1. Let N be an integer > 2, let H be a Hilbert space, and let f G 
H, llv^ll = 1, be chosen such that Sof = f, where {Si)^~^ is a given irreducible 
representation of On on H. Then there is a maximal family {ipn)n>iJ that 
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(3,1) 


ll^nll ^ 1, S^i'n = 0; and 

H = C(^© 

n>l 

Note, every maximal family will satisfy (3.1). 

Before starting the proof, we need a lemma. 

Lemma 3.2. Let-ip G 7-f\{0}. Assumptions as above, including = 0. Then 
the corresponding subspace TL {gp) is invariant under both of the operators Sq 
and Sq. 


PROOF. By (2.6), H {f)) is the closed span of { S^Sjip \ m G Nq, J ^ K [f)) } 
Since Sq%1) = 0, the only non-trivial case to consider is the set of vectors SqSjPj 
when J 7 ^ 0 and J ^ K [f)). But note that SqSjiP = Sj/ip where J' ^ K {fj). 
This concludes the proof. □ 


PROOF OF THEOREM 3.1 Using Zorn’s lemma, we may select a maxi¬ 
mal family {'pn) as stated in the theorem. When it is chosen, our claim is that 
then (3.1) holds. 


We hnish the argument by assuming that the orthocomplement £ of C </? © 

E„>l^(^n), £ = {C99©E„>l^(^n)} = ^ © {C(p©En>l^(^n)}, U 

nonzero, and then derive a contradiction. 

By Lemma 2.2, there are two cases for the projection SqSq\ 

Case 1. S'qS'qI/; ^ Ac, and 

Case 2. S'oS'q|£ = Ac- 

Assuming Case 1, there is some p E C \ {0} such that = 0. Following the 
argument from the proof of Theorem 2.14, we now prove that 

TL (fj) -L TL (pn) for all n > 1; 

and this then contradicts maximality of the family {ipn)- 

Assuming Case 2 , and £ 7 ^ 0, we get 

SiS*p = 0 for alH, 1 < i < iV — 1, and all ip E C. 


This follows from 


N-l 

T. s.s’i’ 


i=0 


Ip. 
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We conclude that S*iIj = 0 for i = 1, 2,..., iV — 1. 

Now use the properties of the vector ip to conclude that 

{Sji/j \ (p) = {i/j \ Sj(p) = 0 for all J G M.. (3.2) 

We have Sop> = Sq(p = (p, and therefore S*(p = 0 for i > 1. As a result, the 
closed subspace spanned by {Sjip | J G WItv} is invariant under On- Since 
9 ? 7 ^ 0, this subspace must be all of H, and we get = 0 from (3.2), which is 
a contradiction. 

This concludes the proof of the theorem. □ 


Proof of Theorem 2.14 concluded. We have chosen to introduce the tools of 
Cuntz algebra representations before the conclusion of the proof of Theorem 

2.14. Once this is available, the reader will easily be able to establish the 
existence of a maximal family of subspaces as in the conclusion of Theorem 

2.14, following the reasoning above (in the proof of Theorem 3.1) using Zorn, 
and in the same way dividing the analysis up into two cases. 


4 A scale-4 Cantor set 


In this section we continue work on harmonic analysis for iterated-function- 
system fractals started in [JoPe98] and continued in [DuJoObb]. Our present 
basis constructions differ from the earlier ones mainly in their emphasis on 
the use of “localized” functions. 


We begin the section with an outline of a new application of our general 
algorithm from Theorem 3.1 to the case of the scale-4 Cantor set X = X 4 on 
the real line M. This Cantor set was studied in [JoPe98]. We begin with a brief 
review of X 4 and its associated measure p = pp. First note that by [Hut81], 
or by direct computation, there is a unique Borel probability measure /i on M 
satisfying 

/ / = ^ (/ / (0 dp ix)+ J f dp (x)^ (4.1) 

for all bounded continuous functions /. Its support X is the unique compact 
subset of M satisfying X = tq (X) U Ti (X) where 


a; 

^0 = 


Ti (x) 


a; -I- 2 
4 


(4.2) 


The set X is sketched in Fig. 2. Specihcally, X = {/^j/4* | /c* G {0,2} }. 
When fractions inside [0,1] are written in base 4, X results as a Cantor set 
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0 


1/4 


1/2 


3/4 


1 



0 1/4 1/2 3/4 1 

Fig. 2. The Cantor set X 4 

from always omitting the two choices 1 and 3. Equivalently (see Fig. 2), by 
repeated subdivision of [0,1], our Cantor set X results when, at each quarter 
subdivision step, we are omitting the second and the fourth open subintervals. 


It is evident that X C [0,1] and that supp (ju) = X where supp (ju) de¬ 
notes the support of p. Furthermore, the measure fi is the restriction of the 

Hausdorff measure Hausi of Hausdorff dimension i. We shall be interested 

2 ^ 

in the algorithms for ONBs in the Hilbert space (p) = {X 4 , ( 14 ). One 

reason for the choice of this particular example is that it appears to be pro¬ 
totypical of IFSs for which (p) has an ONB consisting of complex expo¬ 
nentials ex {x) := Specihcally, the co-authors of [JoPe98] prove that 

{ ca I A G {0,1,4, 5,16,17, 20, 21,... } } is an ONB in (p). (Note that the 
index set for this ONB is the set of all hnite sums 


I /cq -|- kiA / c 24 ^ ks4:^ 


ki e {0,1} I.) 


There are other choices of ONBs, see, e.g., [DuJoOS] and [DuJo06a], but the 
argument for orthogonality is based on a form of factorization of the following 
transform: 


fi (A) := J ex (x) djj, (x) 



The two operators So, Si generating our associated representation of O 2 in 
{X 4 , fj,^) are specihed as follows: 


and 


(Sof) (a:) 


/ (4a;), 

/(4a;-2), 


X ETo (X4), 
X ETi (X4) , 


(4.3) 


(>^ 1 /) (a^) 


/ (4a;), 
-/(4a;-2), 


X ETq (X4) , 
X ETi (X4) . 
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(4.4) 



We have S'oCo = cq, and the formula (4.1) takes the form 


/ fdfx= [ iS*f)ix) dii{x), feCiX). 

J X J X 


We now prove that the representation of O 2 outlined above for the Cantor 
example is irreducible. 

Theorem 4.1. Let Ti = jj,^) be the Hilbert space defined from the 

Cantor measure /i 4 with supp = X 4 . Let Si, i = 0,1, be the operators 

defined in (4.3) and (4.4) above. 

These two operators then generate an irreducible representation of the Cuntz 
C*-algebra O 2 acting on the Hilbert space TL. 


PROOF. We hrst check that the two operators S'j, i = 0,1, are isometries. 
The calculation uses the three facts: 

(i) ro(X)nri(X) = 0, 

(ii) ro(X)Uri(X)=X, 

(hi) p (cF) = y/cpi, {E) for all c G M+ and all E E Bi (X) = the sigma-algebra 

of the Hans 1 -measurable sets. 

2 

Now if / is bounded and measurable, then 

[ d/i = / \f {Ax)\'^ dnix) + [ \f{Ax-2)\‘^ diJ,{x) 

JX J tq{X^ J t\{^X) 

= [ \f{y)f d^M{To{y))+ [ \f{y)\^ dy{n{y)) 

J X J X 

= \ f Iff dy+l- [ Iff dy= [ Iff dy. 

2 Jx 2 JX JX 


The remaining Cuntz property, SqSq -f SiS^ = is immediate from 


y 


- (/i O To"^ +yo rf . 


Set eo = (p = xx = the indicator function of the Cantor set X = X 4 . We 
claim that the vectors Sj(p, J G A4 (= all hnite 0-1 multi-indices), span a 
dense subspace in TL. Since Sq(p = (p and Sip = 0, it follows that p is then 
a cyclic vector for the C> 2 -representation; and that ui^ ■.= {p \ ■ p) \s a. Cuntz 
state. An application of [BJP96, Theorem 3.3] then yields irreducibility. 

For all multi-indices J, K of the same length, we set JK := 'ffiJiki. Recall 
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ji,ki G {0,1} for 02 - An induction now yields 

Xrj(X) = ^^ Y. Sk^. (4.5) 

|Xl=’lJ| 

But by the Cantor construction, we know that these indicator functions 

{ Xtj{x) J G Af I (4.6) 

span a dense subspace in Ti = L^(/i 4 ). The fact that the functions in (4.6) 
are total in is a consequence of the recursive algorithm used in the 

construction of the measure p 4 (see Fig. 2). Moreover, the argument for the 
present special case in fact carries over mutatis mutandis to general IFS con¬ 
structions, regardless of whether the limit measure p is fractal or not; see 
Section 2 above and [HutSl] for further details. 

This concludes the proof of our theorem. □ 

Corollary 4.2. [JoPe98] Consider the Cantor system {X,^) = (X 4 ,p 4 ), the 
Hilbert space Ti = (X, p), and the functions 

Cn (x) := 

n G {0,1,4,5,16,17,20,21,...} 

= { Jo + Ji4 + 32 ^^ + ■ ■ ■ + Jp4^’ I j, G (0,1} } =: A. (4.7) 

These functions form an ONE in TL, and 

77 = C¥P© ^®7f(eJ. (4.8) 

m : 

m odd, 
mGA 

PROOF. The result is a direct corollary of the two theorems. Theorem 3.1 
and Theorem 4.1. In the present application it turns out that K (cm) = 0 for 
all m odd in A. As a result we get that H (e^) is spanned by S^Cm, k G Nq; 
and recall that S^Cm = 6 ^ 4 ^, xn odd in A. □ 
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